On the Beta Function for Anisotropic Current Interactions in 2D 
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By making use of current-algebra Ward identities we study renormalization of general anisotropic 
current-current interactions in 2D. We obtain a set of algebraic conditions that ensure the renormal- 
izability of the theory to all orders. In a certain minimal prescription we compute the /Seta function 
to all orders. 



Left-right current-current interactions in 2D arise in 
many physical systems, for example in the Kosterlitz- 
Thouless transitions and in the study of free electrons 
in random potentials. The chiral Gross-Neveu model 
is the simplest example which is symmetry preserving 
(isotropic) [Q| . Because such interactions are marginal in 
2D, the couplings are dimensionless and there is usually 
no small parameter to expand in order to explore strong 
coupling. 

In this work we determine the /?eta function to all or- 
ders in a certain minimal prescription. This should be 
sufficient for the study of fixed points. The primary tool 
that allows us to isolate the log divergences at arbitrary 
order are the Ward identities for the currents. Kutasov 
performed a similar computation in the simpler isotropic 
case (i.e. non-abelian Thirring model, which is equivalent 
to the chiral Gross-Neveu model) but argued his result 
was the leading order in 1/fc where k is the level of the 
current algebra We beheve our result to be exact. 

We do not expect that there is anything dramatically 
new to learn for the isotropic case. The behavior of mod- 
els like the Gross-Neveu model is very well understood. 
This is in contrast to the anisotropic, i.e. symmetry 
breaking perturbations where one can expect richer phe- 
nomena. In the anisotropic case not all perturbations 
are in fact renormalizable. We find three renormalizabil- 
ity conditions that ensure the theory is renormalizable to 
all orders. Given these conditions are satisfied, we find 
a compact expression for the summation of all orders in 
perturbation theory. 

We know of no reason to expect additional non- 
perturbative corrections to the /3eta function due for in- 
stance to instantons. Our result should be a useful tool 
for exploring strong coupling physics for the many mod- 
els in this class. In this Letter we only report the main 
result and defer applications to future publications. 

Consider a conformal field theory with Lie algebra 
symmetry realized in the standard way |^,^ . It possesses 
left and right conserved currents, J°'{z), J°'{z), where 
z = X + iy, z = X — iy, satisfying the operator product 
expansion (OPE) 
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metric (Killing form) and k is the level. 

We include the case of Lie superalgebras with applica- 
tions to disordered systems in mind. Here each current 
J° has a grade [a] = or 1, and the tensors satisfy: rj'^^ — 

^_^la][b]^ba jabc _ _ |-_ -j [b] [c] yacfc jab _ _ j-_ [a] [fc] j6a 

where f"''"^ — ff'^rf'^ ■ For superalgebras, 77"'' is generally 
not diagonalizable, but we have rjab'ff"^ — (5^. 

The conformal field theory can be perturbed by mar- 
ginal operators which are bilinear in the currents. The 
most general action is 
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where S'cft is the conformal field theory with the current- 
algebra symmetry, and d^^ are certain tensors that define 
the model. The simplest case is that of a single coupling 
where daa = Vaa so that the perturbation is built on the 
Casimir and preserves the symmetry. 

For these models the /3eta function beyond 2 loops is 
prescription dependent. Let /3g = 625^ + 635'^ -I- . . . . The 
prescription dependence is equivalent to a redefinition of 
the coupling g' = g'{g). Let g' = g+cg^ + . . . . One easily 
sees that f3'{g') = ^23'^ + b^g'^ -}-..., so that the one- and 
two-loop coefficients 62 and 63 are universal. Though the 
higher loop contributions are prescription dependent, the 
existence of fixed points is not: if /3g = , then f3'{g') — 
— . For this reason, the particular prescription 
we will adopt is meaningful. 

We will compute the /3eta function in the following way. 
Let X be any arbitrary product of fields and suppose 
we isolate the ultra-violet logarithmic divergences in the 
following way: 



(X) = ^F^(5,loga) / '^{0^{x)X) 



where a is a short distance cutoff (lattice spacing), 
general Fa has an expansion in powers of log a: 



Fa 



(.9) 



n=l 



(5) log 



and similarly for J°-{z). r]°'^ in the above equation is a 



Requiring dFA/d\oga — gives 
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In our models equation (|[) holds due to the Ward 
identities for the currents. This makes the computation 
much simpler that the two-loop computations in per- 
turbed conformal field theory carried out in ||^. The 
order term in the expansion is 

(~)" A A 
j-ffAi ■ ■ • gA^d^lai ■ ■ ■ da^ar. 



X j di...dn{r^r^{i)...r-r-{n)x), (6) 

where J dl = ^ d^xi/2'K. The Ward identities read 

{r^{i)...r-{n)x) 

+ ^{rH2)...r {{)... r"{n)x) , (?) 

Zli 

where Zij ^ Zi — zj, and J means the current is removed. 
(We do not display the factors of (_ in the above 
equation and eq. (||) below.) In order to compute Fa 
one only needs the integrals 
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Stt (p - w) P - z 2w ~ z 



= — log |w — zj , (8a) 
(86) 
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^<5(2)(u;-z), (8c) 



to systematically reduce (^) to one integral. In perform- 
ing this reduction, we throw away two kinds of contri- 
butions: (i) vacuum bubbles proportional to the volume 
J d^x; (ii) logo divergences that factorize into a lower 
order contribution to a*-^^ times a finite part which is 
cancelled by dividing by Gab- The remaining log a di- 
vergences have a rather simple structure. Consider doing 
the integral / dl in (|^) using One is left with 3 kinds 
of terms [fc^], [kf], and [/^] terms: 
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X (^...j r^ii)... ...x'j + z^z^ 

[f] = - /d2 . . . log f-^-^f-^^-^^ 

X {...rr^^{i)...r^j%j)...x) (9) 



Log a divergences arise in the following ways: (i) The [/^] 
term gives a log divergence when the remaining integrals 
give S''^\zij) by repeatedly using (8c). (log|zy [ = loga 
when i = j, (ii) Repeatedly using (86) in the [kf] terms 
can flip 1/z to 1/z or vice versa a number of times. Ul- 
timately, a Ward identity involving another / will give a 
log(zij) like in the term. Then again, the remaining 
integrations must be proportional to (5'-^-'(zjj) to give a 
logo divergence. 

Integrations of log(z.y) with other integrands which are 
not (5-functions give rise to higher powers of log a but are 
not important for the /3-function. We checked explicitly 
at two loops that the log^a divergences are squares of 
one-loop divergences. 

The above contributions can be represented graphi- 
cally, which will be useful for determining their tensor 
structure. Let solid lines denote the J{z) OPEs and 
dashed lines denote the J{z) ones. Denote the /-OPEs 
with arrows. For the holomorphic OPEs for instance: 
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FIG. 1. Diagramatic representation of k- and /-OPEs. 

First consider n = even, which corresponds to an odd 
number of loops. There are two types of contributions: 
Type 2A, 2B. Type 2A is shown graphically in Figure]^: 

n „ m+2m+l 

\ 

m V 

\ 2 3 m-l/ 



FIG. 2. Type 2A diagram. 

The diagrams represent successive integrations from top 
to bottom. The top alternating holomorphic and anti- 
holomorphic lines represent succession of l/z to 1/z flips 
using (86). The bottom alternating lines represent S- 
functions coming from (8c). The rest of the diagram 
represents the /" ''/g'' leading to the log divergence com- 
ing from (8a). The external lines represent the left- 
over J J in (||). In the Type 2A diagrams, m is odd, 
TO = 1, 3, . . . , n — 1, and there are n/2 diagrams of this 
type. The Type 2B diagrams are distinguished by having 
the current that remains after the OPE becoming an 
internal part of the diagram: 




FIG. 3. Type 2B diagram. 



Again to is odd, m = 3, 5, 
diagrams of this type. 



1. There are (n-2)/2 
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At n = odd order, diagrams involve two holomorphic 
/'s or two anti-holomorphic /'s. There are again 2 types, 
lA, IB. The lA diagrams are 
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FIG. 4. Type lA diagram. 

Here m — odd = 3, 5, . . . , n — 2. There are (n — 3)/2 of 
these. A Type IB diagram is distinguished by both /'s 
coming to the same point: 



n— 2 n— 1 




FIG. 5. Type IB diagram. 

There is one diagram of this type. It is the only contribu- 
tion at two loops. For n odd each of the above diagrams 
has a relative" where dashed lines are replaced with 
solid and vice versa. These give the same contribution. 
The total number of diagrams is then (n — 1) for n even 
or odd. 

Each diagram represents a tensor built out of ry"'' and 
which must be contracted with d^^'^- arbitrary 
choice of d^j, does not lead to a renormalizable theory, 
i.e. (@) does not hold. At one loop one needs: 



c 



(10) 



for some structure constants C^^ . This closure of the 
operator algebra of is equivalent to the OPE 
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and it is well-known that C^^ determine the /3eta func- 
tion to one loop [Q. This one-loop condition does not 
ensure renormalizability at higher loops. Systematically 
examining the higher loop corrections we find that one 
also needs the quadratic form d-^ to satisfy 
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The conditions (|10|) and (|13|) guarantee that the theory 
is renormalizable to all orders. The conditions ( p^ ) al- 
ready arise at 2 loops. One can show from the defin- 
ing relations that D,C satisfy D^^ = D^'^,C^^ = 
r<BA riAC nDB _ riAB riDC 



The structure constants D,i? also have a meaning in 
terms of OPEs. Define T'^(z) as 

T^{z)^dtr{z)jHz) (14) 

suitably normal ordered. is a kind of energy-mom- 
entum tensor built out of and the left-moving currents 
only. In the symmetric case is the affine-Sugawara 
stress-tensor up to a normalization. Then one finds 

T^{z)0^{Q) - ^ {2kD^^ + R^D^^) 0^(0) , (15) 

In practice this is the most efficient way to compute the 
RG data. The D term is distinguished from the RD 
term by being proportional to k. As shown below, the 
/3eta function depends only on the combination RD so 
there is no need to compute R separately. The conditions 
under which T'^ actually define a Virasoro algebra ||^ are 
stronger than our renormalizability conditions. 

Finally we can write down the beta function. Let us 
arrange the couplings into a row vector g = (gi, g2, ■ . ■)■ 
It is convenient to define a matrix of couplings D{g), 
^{9)b ~ X^c ^B^ 9c- Then gD^[g) is also a row vector. 
The 4 kinds of contributions (coming from the Type 2A 
through Type IB diagrams respectively) are 
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where D — D(g). (Each diagram comes with a symmetry 
factor of n!/2.) 

Remarkably, the series can be summed in a closed 
form. To display the result in a compact way, for 
any two vectors and v'^ let C(w^,w^) denote a new 
row vector, C{v^,v'^)a = c ^s^c ^A*^- 
g' = .9 (1 - k^D^{g)/Ay\ then 
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-C{g',g'){l + k'Dy4) 



+ '^Cig'D, g'D)D - ^g'DRD , (17) 

where again D ~ D{g). The above /3eta function deter- 
mines how couplings flow with increasing length scale a, 
where large a corresponds to the low-energy limit. 

The simplest possible case corresponds to a single cou- 
pling with dab — Vab- Normalizing the Lie algebra genera- 
tors in the fundamental representation F so that Tr(t°i'') 



cpS'''' = 7]°-'' and using 'q^jfl'^ff = Cadjf^^j, one finds 
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This agrees with previous two-loop calculations for the 
Gross- Neveu model (The usual convention corre- 

sponds to cp — 1/2 with fc = 1.) A three-loop com- 
putation was performed in . 

From the result ([ij) we see that the /?eta function is 
identically zero if Cadj = 0. This occurs for the super- 
algebras osp{2n + 2\2n) and PSL{n\n) 0-|l|]. In the 
work |l^, an all orders /Seta function was proposed in 
the case oi gl{l\l) using a simple scaling argument, and 
the result had a structure of the kind described here. 

The renormalizability conditions ( pl] , p^ are rather re- 
strictive. Consider for instance a completely anisotropic 
perturbation J2a9a.J"'J'^ = Sa5a^° ™ ^ basis where 
j^ab _ gab^ - g ^j^g coupliugs ga,gb are not equal for 
a ^ h. This is not renormalizable in general since the 
operator algebra of does not close: 



C'"(z,z)O''(0) 



(19) 



The algebra only closes if , t'\ is proportional to a 
single generator, since then ff'^f' oc bed and one has 
Cf' = (/c'')^- Adopting a basis for J° based on the root 
system of the Lie algebra, in order to build a renormaliz- 
able theory one needs to include additional interactions 
of the form (a ■ H){a' ■ H) where a, a' are roots and H 
is in the Cartan basis. The resulting theory is renormal- 
izable and one can then compute the RG data C,D,R 
and the f3eta function. Interesting models which are still 
renormalizable can be obtained by equating subsets of 
the couplings consistent with a global sub-symmetry. 

The simplest example is su(2) since it only has one 
Cartan generator. Let us illustrate our main result in 
this example. Let the currents be normalized as 

J3(z)J±(0) ^±-J±(0), J+(^)J-(0) ^ |l + lj3(0), 
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This corresponds to ?f'' = ^S"'', Cadj = 2, and = 1/2. 
Consider the interaction 

J29aO^ = g^{j+J- +J-J+) +g\\J3J3 . (21) 



Using (|Tl|) and (|T|) the RG data are: Dj^-^ = : 
C]^" = = -1, Cf^ = -2, {RD)^^ ^ iRD){ 
(i?J9)j[" — 2. After some algebra one finds 



1/2, 
= 1, 



g± 



P. 



9 J- (.g|| ~ fc.gi/4) 

(l-fc2gi/16)(l + fc5||/4) 
5i(l-A:g||/4)' 
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{l~k^gl/16Y 
When g± = g\\ one recovers 



(22) 



One can use these /3eta functions to study Kosterlitz- 
Thouless flows [|l^ at strong coupling, or the issue of 
symmetry restoration recently addressed in [ p^ . One 
observes flows to the isotropic line g± = g\\ in the regime 
< 9J-,\\ < 4/fc. However there are other features which 
will be reported elsewhere, along with applications of the 
above results to disordered fermions jl^ . 

We hope to publish a longer version of this paper with 
more examples [Q. 
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